Numerical evaluation of a two loop diagram in the cut-off regularization 
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The sunset diagram of A(^* theory is evaluated numerically in cutoff scheme and a nonzero finite 
term (in accordance with dimensional regularization (DR) result) is found in contrast to published 
calculations. This finding dramatically reduces the critical couplings for symmetry breaking in the 
two loop effective potential discussed in our previous work. 
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I. INTRODUCTION 

Dimensional regularization [1] has become an efficient regularization scheme for theoretical calculations in high 
energy physics. It makes analytical results more available, and its gauge invariance almost decisively defeated the 
other schemes (say cutoff scheme) in QCD and other gauge theories. 

However, things are not so simple. It is very natural to anticipate that DR might fail sooner or later somewhere 
as it is just a special regularization but never a true theory of the short distance physics. Such event is now known: 
the MS scheme does not provide a physical prediction for low energy nucleons scattering in case of long scattering 
lengths [2] in nonperturbative context. To yield meaningful results, unconventional subtraction is devised and now 
known as PDS [2]^. Since this is obtained with hindsight, we should not discard the cutoff scheme without further 
exploiting its virtues. There have recently been some further investigations over the relations between the divergences 
in DR and cutoff schemes and the associated limitations [4,5]. The conventional virtues of the DR scheme were shown 
to be due to an implicit subtraction [4]. In nonperturbative contexts involving multiloop corrections the DR might 
even make renormalization infeasible [5] . 

Traditionally, it is very hard to obtain analytical results for multiloop amplitudes within the cutoff scheme. In this 
report we show that one could achieve this goal numerically at a pretty good precision. We will make use of the 
dynamical symmetry breaking of discussed in [6] to illustrate the numerical evaluation of the finite constant term 
in the two loop sunset diagram. Thus the report is organized as follows. First we quote some results of Ref. [6] in Sec. 
II as preparation. In Sec. Ill we present our strategy and procedures for doing the numerical analysis of the problem 
and its possible prospective significance. The result is also given there. In Sec. IV we revisit the problem introduced 
in Sec. II to show the significance of the new result obtained in Sec. III. The summary will also be given there. 

II. RENORMALIZATION OF THE TWO LOOP EFFECTIVE POTENTIAL 

In [6] , the dynamical symmetry breaking of massless Xcj)^ model with Z2 symmetry was studied in two loop effective 
potential. The algorithm in use is well known according to Jackiw [7] , 



s M' + -J„(!i) + 3A/ji(!i) - 48A»0V2(n), (1) 



' ^ J (27r)8 (fc2 + n^){p + n^){{k + ly^ + n^) ' 

Here we have Wick rotated all the loop integrals into Euclidean space. The loop integrals and the effective potential 
had been calculated in literature [7,8]. For comparison, we list out below only the results for the sunset diagram 
defined by /2(f^) in Eq.( 3) , the only two loop diagram with overlapping divergences. 
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There are other schemes that could yield similar predictions as PDS did, see, e.g., Ref. [3]. 
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A. Loop amplitudes in DR and cutoff schemes 



From [8] we find that 

M^^/f = + (3 - 2L)e + {2L' -6L + 7 + 6S- ^C(2))e^} (4) 

with S = J2n=o (2+3n)^ '^^ L + 7 - In 47r and L = In ^. While from [7], 

Note that the finite double logarithmic term (leading term) disagree between the two regularization schemes, the 
disagreement is removed only after all sub divergences are removed [5] . 

We note in particular that the term proportional to il^ in the two loop integral is not explicitly given in [7] , and it 
is this finite term and its numerical estimation that is our main concern in this paper. 



B. Prescription dependence 



'{21) 

Weinberg (CW) [9] and IJ,\ (defined in [6], Appendix B) for V^^^{Cl) with the renormalized potential taking the 



The renormalization were performed in Ref. [6] in several prescriptions, MS for Vf2n^(^) [^]' J^ckiw [7], Coleman- 

/einberg (CW 
following form 



where L = L, In In and In jttt^ ir respectively MS, u\, Jackiw and Coleman- Weinberg schemes. 

In all the above formulas the scheme dependence of field strength and coupling constant are understood. The explicit 
intermediate renormalization prescriptions dependence of the effective potential expressed by a is listed in Table I. 
The dynamical symmetry breaking solution in the two loop effective potential defined in Eq. (6) could be obtained 

from the first order condition ^^^^'^ P^*^^^ = and the solutions read [6], 



1 

(t>\{\; [m,q;]) = Y^exp j - 
whose existence requires that 



47r2 / 1 47r2 



(7) 



1 , , 477^ 

That means not all schemes arc consistent with symmetry breaking. In addition, for the symmetry breaking solutions 
to be stable, the coupling must be further constrained [6], 

47r2 

' ^ " V4- 36.- 27- 1 ' ^> 

Both Eq. (8) and (9) are summarized in Table II. 

If the finite term proportional to f2^ in l2{^) existed, the a in fjb\ and Jackiw prescriptions would consequently be 
replaced by — 2 — ^ and — | — ^ in the following parameterization of Eq. (5)^, 



^In Coleman- Weinberg prescription the definition of a = ^ is not altered due to its special renormalization condition [7,9]. 
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It is not difficult to sec tliat a different value of a means different symmetry breaking status, hence different physics, 
a nonperturbative scheme dependence problem as emphasized in ref. [6]. It is therefore important to determine the 
value of C in cutoff scheme in order to be sure that the Jackiw and ^ prescription are really consistent symmetry 
breaking in two loop effective potential. Thus the rest of this report will be devoted to estimate this constant term 
numerically. 



III. STRATEGY FOR NUMERICAL ESTIMATION OF C 

First let us put the sunset diagram into the following form after integrating out the angular variables, 

with M = fl^,x = k'^,y = and /? = A^. The asymptotic form of this integral as cutoff tends to infinity reads 

HM) = |2/3 -^^J + 3Mln j + Cm| + o(M//3) 

^7^^)(M) + ^+o(M//3). (12) 

In order to determine the constant C we compare I2 (M) and 1^°''^'^^ (M) after we evaluated the former numerically. 

The feasibility of numerical treatment of an apparently divergent integral is guaranteed by the observation that the 
cutoff need not be very large. Eq.( 12) tells us that the analytical source of error comes from The magnitude 
of the constant C should be of order 1 or 10 in the limit 00. Thus taking to be 10^'* or 10^^ should be 

sufficient for our purpose. There is a main obstacle in the practice for numerical estimation of divergent integral 
in the cutoff scheme: the presence of power law pieces, which consume most of the computer's capacity in order to 
get a finite number at a precision of 10""*. The larger the cutoff is chosen, the less precision capacity is left over 
for finite numbers. We must get rid of such pieces. The remaining logarithmic terms save most of the computer 
capacity for the finite numbers, as even taking the cutoff as large as 10^°M will only yield a double log term at order 
~ 100 In 10 ^ 230, which consumes almost no precision capacity. This goal could be readily achieved at least in two 
ways: for the sunset diagram, the power law term ~ 2/3 could be removed either by subtracting it explicitly in the 
integrand or by differentiating the whole integral with respect to M. In fact we will report both approaches' result. 

There is another subtle point that need careful consideration, the step length for the numerical integration. Taking 
the step too small will improve the precision analytically but cost a lot of time and might accumulate numerical 
errors, especially in personal computer. Usually a divergent integral has a smooth integrand except in a finite range 
where the integrand is steeply changing. The step must be small enough within such steep regions. The last concern 
is to spare both the computer time and numerical error we should perform the analytical calculation as far as we can. 
Thus we will first perform one integral analytically (this is feasible in most cases) while leave the remaining one to 
computer. 

In the following we will describe one approach of our numerical treatment in detail, the one with differentiation 
with respect to M. That is, to extract C we just consider the following double integral which at most possesses 
logarithmic divergence, 



dM (47r)4 \ ' " ' 2' ^ ^'j dM (47r)4 

dx /-^ dy 

lo {x + Mf Jo V+M' 



h^l TTTTJ^l ^i^ + y + M), (14) 



'^-L ¥TWL ^Vi^^y^Mr-.xy, (15) 



lo x + M Jo y + M 
J _ dx f" dy x + y + M 

^ Jo x + M Jo y + M ^(^x + y + MY - Axy ' 
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It is not difficult to sec that the 2/? piece is removed in all the four integrals Ii, h, h and there is at most double 
log terms. Then it is also easy to perform two folds integrals once. Moreover, Ii and I3 can be done analytically at 
all. Then the results read 

^JTM~^rM{M + p) ' ^^^^ 
X . I a — xY . I Y 



h 


, 13 + M ( P + M 
= in — T7 — i in — — — 
M \ M 


h 


_ dx 1 




~ Jo x + M] 




h 


2I3 + M 
= M ' 




h 


r'^ dx J 


— ^sinh" 


~ Jo X + M] 


2xY 


+ y2 a = x^ 


+ 4:xA4. 



^ ^^^^ + sinh- ^ \ 11^1^+^ , (19) 



(20) 

-sinh-^|||-=r^ (21) 



with R = a — 2xY + Y'^, a = x^ + AxM. There is a seemingly singularity in the integrands (~ as a; ^ 0) in 

both I2 and J4 which does not materialize after one completes the operations on the variable y, 

' f "^lirt:^! ,^ ^^u+V^l ^^s + ./srTl\ 

x + M\ y 1'^-^ Va: + 4M d + VdPTl t + V^F+l \ 



-L 



dx , M + Vu^+T s + Vs^ + 1 I 



x(a; + 3M-/?) + 3xM fi + M-x M-x 

— ; , V = « = + = 

2{f3 + M)y/^ 



(24) 



The only necessary numerical works to be done are the two one fold integrals defined in Eq.( 22, 23). As the integrand 
become steeper as x is closer to the origin, thus to reduce the error from numerical integration we separate the integral 
into several intervals in order to match between the step size and the steepness of the integrand. Taking M to be 1 will 
make the numerical integration easy and the C has been estimated with two cutoff scale, /3/M = A?/Q? = 10^, 10^. 
The match between the interval and step size is shown in Table III and IV respectively for f3/M = 10* and P/M = 10^. 

Subtracting — ^-gfj' — from the numerically obtained integral defined in Eq.( 13) we find the constant C: C = 
4.13048 for /?/M = lO'^ and C = 4.15412 for (3/M = 10^. The interval (0,10"'^] is not included as its contribution 
could at most be of order 10~*. Noting that one cutoff is order greater than the other, the agreement of the two cases is 
striking. Moreover, if one try another approach, i.e., by directly subtracting /2's integrand a term that upon integration 
will yield the 2/3 piece, one could also obtain an estimate of C with the results read: for (}/M = 10'', C = 4.14128 
with a homogeneous step 10~*, C = 4.14274 with a homogeneous step 10~^, C = 4.156178 with a homogeneous step 
10-2; for P/M = 10^, C = 4.160739 with a homogeneous step 10"^ C = 4.1633 with a homogeneous step 0.005; for 
P/M = 10^, C = 4.162993 with a homogeneous step 10"^. Combining the two approaches' results we could safely 
conclude that C is not zero. To be more conservative we anticipate that the true value should lie in the interval 
[4,4.3]. 

Before ending this section, we should point out and stress that in order to evaluate divergent multiloop integrals 
numerically in cutoff regularization: (i) one does not need a very large cutoff scale, (ii) the power law term can 
be easily removed in order to save capacity for more important task (like determination of the finite parts). The 
numerical workload could further be reduced if one varies the step size with steepness of the integrand. 



IV. NEW PREDICTIONS AND SUMMARY 



With the new constant we could reevaluate the critical couplings for symmetry breaking with the result summarized 
in Table V and VI. We take C to be 4 when calculating the critical couplings. It is not difficult to see that the 
critical couplings in /i^ and Jackiw prescriptions are dramatically lowered, or the symmetry breaking might take 
place at much lower values of the coupling. The interesting thing is, the MS critical coupling, which is smaller than 
the ones in the prescriptions based on old calculation of the sunset diagram [7] , now becomes a larger one than the 
obtained with the new numerically determined constant C, and the larger the C is, the smaller the critical couplings 
are. In this sense, the cutoff regularization is preferable to DR, just like the case in the EFT applications in nucleon 
interactions [10]. 
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Generally one would expect the critical coupling to be close to 1 instead of being very large (like 10.7) for symmetry 
breaking to take place. Therefore it is a pleasing finding that constant C is larger than and is no less than 4. If 
C = 20, a reasonable value, the critical coupling could be 1.5238 and 1.5542. Of course, as the potential is only 
calculated at two loop level, the precise value of the critical coupling could not be taken very serious. However, one 
should still be pleased to see that the critical couplings are dramatically lowered after more careful evaluation in the 
regularization scheme that is widely held as inferior to DR, even though the computation is not exact one. 

In summary, we recvahiated the sunset diagram in A(/)^ in the cutoff regularization half analytically and half 
numerically and found that the finite local constant, which is usually taken to be zero, is not zero and this constant 
could dramatically lower the critical couplings for dynamical symmetry breaking in the two loop effective potential 
in a number of renormalization prescriptions. The important byproduct is that the numerical calculation could be 
efficiently done in cutoff regularization for divergent multiloop diagrams and we have illustrated that the cutoff scale 
need not be too large. 



TABLE 1. a in various schemes 



Scheme 


Q 


MS 


-2.6878 




-2 


Jackiw 


5 
4 


Coleman- Weinberg 


49 
3 



TABLE 11. Critical values of A 



Scheme 


Acr 




MS 


4.368 


5.2024 


lA 


5.1152 


6.5797 


Jackiw 


6.5797 


10.698 



TABLE 111. Match 1 for (i/M = 10* 



Interval 




[lO"^lO"^] 


[10-\10"] 


[10", 10^] 


[10\10^] 


[lo^lo^] 


[lo^lo*] 


Step Size 




10-' 


10"" 


10"" 


10"" 


lO"-" 


lO"'' 



TABLE IV. Match 11 for /3/M = 10^ 



Interval 


[10-®, 10"''] 


[lo-^lo-^] 


[io-\io°i 


[10°, 10^] 


[10\lO^] 


[lo^lo^] 


[lo^lo*] 


[10*, 10^] 


Step Size 


10"° 


10-' 


10"° 


w-° 


10-* 


10-" 


10-" 


10-" 
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TABLE V. New a's 



Scheme 


(I 


MS 

lA 

Jackiw 
Coleman- Weinberg 


-2.6878 
-7i 

16| 



TABLE VL New critical values of A 



Scheme 




Acr 


MS 


4.368 


5.2024 


A 


2.5684 


2.7181 


Jac;kiw 


2.71cSl 


2.CS967 
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